Unitary polarized Fermi gas under adiabatic rotation 
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We discuss the effect of an adiabatic rotation on the phase separation between the superfiuid 
and normal component of a trapped polarized Fermi gas at unitarity and zero temperature, under 
the assumption that quantized vortices are not formed. We show that the Chandrasekhar-Clogston 
limit ni/n-\ characterizing the local polarization in the normal phase at the interface is enhanced 
by the rotation as a consequence of the centrifugal effect. The density profiles (local and column 
integral) of the two spin species are calculated as a function of the angular velocity for different 
values of the polarization. The critical value of the angular velocity at which the superfiuid exhibits 
a spontaneous quadrupole deformation is also calculated for the unpolarized case. 

PACS numbers: 03.75.Ss, 05.30.Fk, 47.37.+q, 67.90.+Z 



I. INTRODUCTION 

Recent experiments 0, [H, S 0, H| have shown that a 
polarized Fermi gas at unitarity and zero temperature 
undergoes a phase separation between a central unpo- 
larized superfiuid component and an external polarized 
normal gas. Experiments, where surface tension effects 
are not important and the local density approximation 
(LDA) is applicable, have revealed the occurrence of a 
critical value of the polarization 



P = 



(1) 




of the gas above which the superfiuid component disap- 
pears. It is found that such a critical value is P c ~ 77%. 

The study of spin-polarized Fermi gases has been the 
object of several theoretical papers. In [|| Q the equa- 
tions of state of uniform matter for the superfiuid and 
for the polarized normal phases, calculated with ab ini- 
tio Monte Carlo simulations, have been employed within 
a local density approximation to treat the effect of the 
harmonic trapping in the unitary regime [H . These calcu- 
lations were proven to be very efficient not only in repro- 
ducing the experimental value of the critical polarization, 
but also the density profiles of the two separate spin com- 
ponents. Therefore they provide an accurate and consis- 
tent description of the phase separation exhibited by the 
unitary Fermi gas at zero temperature. In particular, 
the discontinuity characterizing the spin-down density at 
the interface between the superfiuid and normal compo- 
nents, as well as the typical knee revealed by the column 
density of the same spin-down component, are dramatic 
features reproduced with high accuracy by theory. These 
calculations have pointed out the crucial role played by 
the interactions in the normal phase Q- 

The purpose of this paper is to explore the rotational 
properties of these polarized Fermi gases (see the sketch 
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FIG. 1: (Color online) The typical shell structure of the 
trapped system consisting of a superfiuid core (red) sur- 
rounded by a partially (green) and fully polarized (blue) nor- 
mal shell for a polarization P = 44% and f2 = 0.5lo±. The 
superfiuid is squeezed in the radial direction while the normal 
part exhibits the bulge effect of the rotation. 



in Fig. []}. It is well known that the response to a trans- 
verse probe, like the rotation, is a crucial tool to test 
the superfluidity of a system. In a trap, this can be 
achieved by rotating the confining potential. While a 
normal gas rotates in a classical rigid way, a superfiuid 
features a different behavior. The phenomena exhibited 
by superfluids include the quenching of the moment of 
inertia at small angular velocities and the appearance 
of quantized vortices at higher velocities. Vortices have 
already been observed in these polarized configurations 
and shown to disappear for high polarizations. In a re- 
cent paper Q we have predicted that a trapped rotating 
Fermi gas at unitarity exhibits a further interesting phe- 
nomenon, associated with the breaking of superfluidity 
in the external region. Indeed, atoms prefer to be in 
the normal phase because of the energy gain associated 
with the rotation. This mechanism of the depletion of 
the superfiuid due to the rotation has also been recently 
confirmed within BCS mean-field theory The oc- 

currence of this phenomenon requires proper conditions 



2 



of adiabaticity in the ramping of the rotation of the trap 
in order to avoid the formation of vortices, a condition 
that has been already successfully realized in rotating 
Bose-Einstein condensates [ll| . 

In a polarized Fermi gas the phase separation be- 
tween the superffuid and the normal component is al- 
ready present in the absence of rotation, so that the ef- 
fect of the rotation is expected to further enhance the 
mechanism of depletion of the supcrfluid with new visi- 
ble consequences on the density profiles of the two spin 
components. The explicit investigation of these features 
represents the main goal of the present paper. 

This paper is organized as follows. In Sec. |TT] we in- 
troduce our model and investigate the effect of rotation 
on the phase separation between normal and supcrfluid 
within the LDA. The results are discussed in Sec. IIIII In 
Sec. IIVI we calculate the critical frequency for the break- 
ing of axi-symmetry of the supcrfluid core and finally in 
Sec. |V]wc draw our conclusions. 



II. POLARIZED FERMI GAS IN A ROTATING 
TRAP 

As already mentioned throughout the paper we assume 
that only two homogeneous phases are possible: A po- 
larized normal phase and an unpolarized superfluid one 

a 

The normal state of a polarized Fermi gas at unitarity 
has been first introduced and discussed in [1, [HI, [H| . Fig- 
uratively it can be understood of spin-f atoms 
to which constantly spin-| atoms are added. Since at 
unitarity there are no interaction length scales, the en- 
ergy of the system at zero temperature can be written in 
terms of the Fermi momentum /cft of the majority com- 
ponent and the concentration x = n j /n j of the minority 
atoms. As soon as the concentration reaches the critical 
value x c = 0.44 the system starts to nucleate a supcrfluid 
region. 

In the unitary limit of infinite scattering length the 
dependence of the energy of the normal phase on the 
concentration x can be written as @, 0] 



E N (x) 



= -E FX I 1-Ax 



? 5/3 + Bx 2 
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E F ^g{x) = e N (x), 



(2) 



where iVt is the total number of j atoms, Ep^ = 
?j- 2 /2m(67r 2 n|) 2 / 3 is the ideal gas Fermi energy, and m 
the atomic mass. The values of the spin-down single- 
particle energy A and of the effective mass to* have been 
calculated both analytically and numerically [H, HH, [l(| , 
the parameter B is known by fitting Monte Carlo re- 
sults. The most recent Monte Carlo calculations give 
A = 0.99(1), to* = 1.09(2), and B = 0.14 [13. The 
paramctrization © reproduces the Monte Carlo results 
for the energy of the normal state not only in the low 



x regime, but also for large values of the concentration 
parameter. 

On the other hand the equation of state for the super- 
fluid phase is simply given by 



Es 
N s 



?S^(67r 2 ns) 2 / 3 ^6 S (n s ), 
5 zm 



(3) 



where As is the number of atoms in the supcrfluid phase, 
ns the supcrfluid density (equal to the spin-up and spin- 
down density) and the interaction parameter £s = 0.42 is 
known from ab initio quantum Monte Carlo simulations 

MM- 

We consider a polarized Fermi gas at unitarity confined 
by a harmonic potential V(r) = m(uj 2 x 2 +uj 2 y 2 +uj 2 z 2 )/2, 
rotating with angular velocity O along the z axis. We 
study the problem in the rotating frame of the trap, 
where the potential is static and the Hamiltonian con- 
tains the additional term — HLz- We stress once again 
that the response of a superfluid to an external rotation 
is multifaceted and can depend on whether one ramps up 
the angular velocity very fast or in an adiabatic way. In 
the following we will not take into account the formation 
of vortices. This is best ensured by assuming that the 
angular velocity is ramped up adiabatically. 

In the LDA the grand canonical energy of the rotating 
configuration at zero temperature takes the form 



dr e(n T (r),n x (r)) + V(r) 
+ —mv 2 — mQ(r x v) z j n(r) 



- / dr[ / u°n T (r)+/i°n i (r)], 



(4) 



where e(nj(r), rij(r)) is the energy density per particle 
depending on the n-|-,i( r ) densities of the two spin species, 
v is the velocity field, fi® and /Xj 1 are the chemical poten- 
tials of the t and J, particles, and n(r) — ri|(r) 4-n^(r) is 
the total density. 

We assume that the phase separation in the trap man- 
ifests as the formation of an inner unpolarized superfluid 
core occupying the region r < Rs(0, 4>) surrounded by an 
external normal shell, which is confined to Rs(9, <j>) < r < 
Rn{0,4>)- Here, we term Rs(9,(f>) the interface separat- 
ing the superfluid from the normal phase and Rn{0,4>) 
the Thomas-Fermi radius of the gas where the density 
vanishes. Thus, the integral ([4]) splits into two parts 



E 



2 / dr 

'r<R s 



( es(n S (r)) 
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V(r) + -mvl - mn(r x v s ) z ) n s (r) 



dr 



Rs<r<B. N 

M°n T (r) - (i,°ni(r) 



e N (z(r))n T (r) 



3 



V(r) 



mfi(r x v N )z ) n(r), 



(5) 

where /ig = +/iJ)/2 is the superfluid chemical poten- 
tial and n^i(r) the | and J, densities in the normal phase. 
In the above equation we have distinguished between the 
velocity fields v§ and in the superfluid and normal 
phases, respectively. 

To find the equilibrium conditions, we minimize the en- 
ergy with respect to the densities, to the velocity fields as 
well as with respect to the border surface. In the case of 
the superfluid the velocity field obeys the irrotationality 
constraint and can thus be written as vs = V<F Varia- 
tion of the energy with respect to the velocity potential 
<i> yields the continuity equation 



V • n s (V$ -fixr)=0, 



(6) 



while the variation with respect to the superfluid density 
ns yields the LDA relationship 



^£s^(6^ 2 n s ) 2 / 3 + U s (r), 



(7) 



where U s (r) = V(r) + ^mv 2 — mfl(rxvs)z is the effective 
harmonic potential felt by the superfluid. 

Using the same procedure for the normal part (with- 
out the irrotationality constraint) we get vn = f! x r, 
i.e. it rotates rigidly. The variation with respect to the 
densities gives the LDA expressions 

M? = (g(x) - £(6-V) 2/3 + V* ('), (8) 



4 



y {x )^(6^n^ + V^r), (9) 
5 2m 



where the effective potential U N (r) felt by the particles 
in the normal phase is now squeezed due to the rigid 
rotation according to (uj^) 2 = w 2 - ft 2 , (uj^) 2 = u\ - fl 2 . 

By varying the energy ([5]) with respect to Rs{0, </>) we 
eventually find the equilibrium condition for the coexis- 
tence of the two phases in the trap. This is equivalent to 
implying that the pressure of the two phases be the same 



2 de s 
n s 



'dm 



S / r=R s 



2 <9e N (» 



■ n\ny 



dni 



t I / n> 

I / r=R s 

(10) 

Using the expressions for the energy densities ([2]) and (J3]) 
we obtain an equation for the density discontinuity in the 
trap given by 



2& 



n s (R s ) ~ \g(x(Rs)) 



3/5 



= j(x(R s )), (11) 



where x(Rs) is the local concentration at the interface. 
Combining the above equation with Eqs.([7]) and ([8]) we 
find the useful relationship 

(4 - U s (i? s )) = j(x(R s )) (//? - U N (i? s )) , (12) 



determining the surface Rs(9,(j)) separating the super- 
fluid and the normal part. 

Eventually using also Eq.Q we find an expression 
which implicitly defines the concentration as a function 
of the position at the interface 

g(x(R s )) + hl-x(R s )]g'(x(R s )) 
5 

- (2^ s ) 3 / 5 [g(x(R s ))] 2 ^ 
_ 2(V S (R S )-V N (R S )) 



Ef\(Rs) 



(13) 



In absence of rotation V s = U N = V the solution of 
Eq.® yields the value x(R s ) = 0.44 0. This value 
coincides with the maximum concentration achievable in 
the normal phase of uniform matter before phase separa- 
tion. Since the rotation affects differently the potentials 
V" and y N , the value of the concentration depends now 
on the angular position of the interface. It ranges from a 
minimum value x[i?s(0i 0)] — x c = 0.44 along the z axis 
(where the effect of the rotation is vanishing) to a max- 
imum, fl dependent value x[Rs(n/2,<j))] in the xy plane 
(where the effect of the rotation is largest). The effect of 
the rotation is thus to enhance the average value of the 
concentration in the normal phase and hence to favour 
the depletion of the superfluid (sec Fig. [TUland discussion 
below) . This is physically understood by noticing that in 
the rotating frame the atoms in the normal part gain the 
energy |mv^(r) due to the centrifugal force. Thus, the 
energy of the normal part can become smaller than the 
value in the superfluid Q. The main effect is to change 
the critical concentration x c at the interface (see Fig. [2] 
below) . 

Notice, however, that the critical global concentration 
P c for the system to start nucleating the central super- 
fluid core is not affected by the rotation and keeps the 
nonrotating value P c = 77% (see also Fig. [3]). The reason 
for that is easily explained within the local density ap- 
proximation used here. Just above P c the system is com- 
pletely normal and the only effect of the rotation is the 
squeezing of the transverse trapping frequencies, while in 
the z direction the system remains unaffected (see Eqs.((51 
|9])). Since at P c the superfluid is nucleated at the center 
of the trap where centrifugal effects are absent, the lo- 
cal condition for equilibrium between the superfluid and 
the normal component is the same as without rotation. 
The calculation of the critical polarization proceeds then 
in exactly the same way as without rotation, but for a 
simple rcscaling of the trapping frequencies which has no 
effect on the value of P c . 



III. RESULTS 

In the following we will assume lu x = lu v = lo±_ and we 
consider the solution vs = 0, thus a nonrotating axisym- 
metric superfluid and consequently V s = V. In this case 
the local concentration depends only on the polar angle 



4 



and the supcrfluid radius takes the form 
2 (Mg- 7 (x(Jfe))M?) 



m (l-j{x(R s ))) 
( w 2 cos 2 + w 2 sin 2 6> 

7(s(fls)) 
l- 7 (a;(i?s)) 



SI 2 sin 2 i 



(14) 



while the Thomas-Fermi radius i?N of the normal gas is 
fixed by the condition ^ = V N (r) yielding 

ij2 = _n ( w 2 CQS 2 g + ^2 gin 2 _ ^2 sin 2 #) " 1 > i? 2 (6>) . 

(15) 

The values of the chemical potentials are fixed by the 
normalization 



/ dr ras(r) + / 

Jr<R s Jh 



dr n T (r) = N h (16) 



fl s <r<fl N 



and 



[ dr n s (r)+ f dr n^r) = N±. (17) 

Jr<R s JR s <r<R N 
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FIG. 2: (Color online) Concentration ni/nj for P = 44%, 
SI = (black small dashed), SI = 0.3tJx (orange dashed), and 
SI = 0.5cj± (turquoise solid) as a function of the polar angle 



In Fig. [5] we plot the concentration riy / nj at the in- 
terface for P = 44% as a function of the polar angle 8 
for different values of the angular velocity. The figure 
clearly points out the increase of the concentration when 
one moves from the z axis to the xy plane. 

Complementary to this is Fig. [3J where we show 
the concentration ni/n-\ for a fixed angular velocity 
SI = 0.5wj_ and different values of the polarization as 
a function of the angular position. This clearly reveals 
the nature of the normal state by highlighting the two 
extreme and singular cases P = 77% and P = 0%. As 
already pointed out, at the threshold value P = 77% for 
the nucleation of the superfluid, the rotation does not af- 
fect the value of evidenced by the solid black line in 
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FIG. 3: (Color online) Concentration nj/n-f at the superfluid- 
normal interface as a function of the polar angle 6 for Q. — 
0.5uj± and different values of the polarization: P = 77% 
(black solid), P = 44% (turquoise dashed), P = 10% (pink 
dotted-dashed), P = 2% (green long dashed) and P = 
(black thin dashed). 



Fig. On the other hand, in the case that P = 0% the 
critical concentration is constant and singular n^/n-f = 1 
for all angles but 8 = 0. For all other values < P < P c , 
the rotation has a considerable impact on the local value 
nj./n T . 




FIG. 4: (Color online) Radii in units of i?| versus angular 
velocity for 9 = tt/2 of the superfluid (red), n-\ (blue), and ny 
(green) for a polarization P = 44%. 

The density profiles exhibit a typical shell structure. 
In Fig. [J] we plot the radii of the superfluid (red), | 
(blue), and J. (green) component in units of the Thomas- 
Fermi radius of an ideal gas versus the angular veloc- 
ity Vl/uj± for a polarization P = 44%. While the super- 
fluid radius decreases until the supcrfluid core completely 
vanishes at SI = uj±_ , the Thomas-Fermi radii of the j and 
I component diverge for SI = lu±_ as a consequence of the 
centrifugal effect. It is curious to see that while the nor- 
mal part exhibits the typical bulge effect, the supcrfluid 
behaves in the opposite way. In fact, its radial size be- 
comes smaller than the axial one as a consequence of the 
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depletion caused by the rotation, with consequently in- 
version of the behavior of the aspect ratio R±/Rz (see 
Fig.©. 
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FIG. 5: (Color online) Column densities of the raj (blue), 
and nj. (green) component in a spherical harmonic trap for 
a polarization P — 44% and Q = (dashed lines) and Q = 
0.5lu± (solid lines). 
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FIG. 6: Difference of the column densities of the n-f and n± 
components in a spherical harmonic trap for a polarization 
P = 44% and O = (dashed line) and O = 0.5lo± (solid line). 

It is worth mentioning that at large enough angular 
velocities the system exhibits solutions which break the 
axial symmetry [2fj| . Such critical value is predicted to 
be fi cr ~ 0.5w_i_ as we discuss in detail in Sec. IIV1 The 
results in Figs [U GIB and QT] for fl > 0.5uj± (dashed ver- 
tical line) correspond to the axial symmetric solution of 
the problem. 

In an experiment the effect of phase separation as well 
as the radius of the supcrfluid are best revealed as a knee 
in the in situ column density n CTj 2D(p) = / dz n a (r), 
with a =|, J,. These observables can nowadays be mea- 
sured with high precision using phase-contrast image 
techniques. We expect that the position of the knee for 
a fixed polarization will depend on the angular velocity. 
This is clearly shown in the column density of the ma- 
jority (j, blue) and the minority (J., green) components, 



Fig. El for fi = (dashed) and £1 = Q.5uj± (solid), as 
well as in the density difference, Fig. [51 The knee is a 
direct consequence of the discontinuity exhibited by the 
three dimensional density shown in Figs [7] and [HI where 
we plot ng (red), n-j- (blue), and nj, (green) in a spherical 
trap for 6 — and 9 = tt/2, respectively. The densities 
and the radial coordinate have been renormalized with 
respect to the central value of n® and the Thomas- Fermi 
radius i?° of an ideal gas. In accordance to the results for 
the superfluid radius shown in Fig. 3] the discontinuity 
in the density takes place at a smaller value of the radius 
compared to the nonrotating configuration. It is worth 
noticing that this knee is also exhibited by the total den- 
sity of the unpolarizcd gas (Fig. [5]), reflecting the density 
discontinuity produced by the rotation Q . 
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FIG. 7: (Color online) Density profiles for 6 = of the super- 
fluid (red), n-f (blue), and nj (green) in a spherical harmonic 
trap for a polarization P = 44% and ii = 0.5ojx in units of 
the central density of the noninteracting gas (dashed line). 
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FIG. 8: (Color online) Density profiles for 9 — tv/2 of the 
superfluid (red), n-f (blue), and ni (green) in a spherical har- 
monic trap for a polarization P = 44% and ii = 0.5cu± in 
units of the central density of the noninteracting gas (dashed 
line) . 

For 8 = the densities and jump from the 
supcrfluid value ng to the values nj ~ 1.01ns and 



6 



ni = x c n^ ~ 0.44ns as one enters the normal phase, pre- 
cisely as in the nonrotating case Yet for 9 = tt/2 the 
behavior is different and, in particular, for = 0.5uj± the 
densities jump from ns to nj ~ 0.99ns and nj_ ~ 0.55ns- 
The smaller relative jump with respect to the nonrotating 
case reflects the smaller polarization (higher concentra- 
tion x) exhibited by the rotating normal phase at the 
interface. 
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FIG. 9: (Color online) Total column density of the unpolar- 
ized system in a spherical harmonic trap for Q = (dashed 
red line, superfluid at rest). For Q — 0.5a; x the system con- 
sists of a superfluid core surrounded by a normal shell where 
n T = ni. 

Further insight on the effect of the rotation is provided 
by the depletion of the superfluid. In Fig. [TOJwe plot the 
ratio between the number of particles in the superfluid 
Ns and the total number TV (superfluid fraction) as a 
function of the angular velocity. This effect is especially 
pronounced for small polarizations (Fig. [TO] black line, 
P = 0) since the depletion for higher polarization is large 
already in the nonrotating case. 




FIG. 10: (Color online) Evolution of the superfluid particle 
number for different polarizations P = 0% (black solid ), P = 
10% (pink dot-dashed) , P = 25% (orange large dashed) and 
P = 44% (turquoise small dashed) as a function of tt/u}±. 

Finally in Fig. [TT] we plot the angular momentum 




FIG. 11: (Color online) Angular momentum Lz in units of 
the rigid value for different polarizations (black solid P = 0, 
pink dot-dashed P = 10%, orange large dashed P = 25% and 
turquoise small dashed P = 44%) as a function of fi/w_i_. 



Lz = Q j dr(x 2 + y 2 )n^(r) of the system for different 
polarizations. For an axisymmetric configuration the su- 
perfluid does not carry angular momentum which is then 
provided only by the normal component. Hence, the 
more particles are in the normal part, the stronger the 
response of the system to the rotation. 



IV. QUADRUPOLE INSTABILITY INDUCED 
BY THE ROTATION 

In the preceding section we have shown that the ef- 
fect of the rotation applied to a polarized Fermi gas at 
unitarity is to enhance the depletion of the superfluid. 

A further important effect is the spontaneous breaking 
of axi-symmetry. In fact, one expects that above a certain 
critical value £l cr the system exhibits a surface energetic 
instability, undergoing a continuous shape deformation, 
similarl y to what happens in Bose-Einstein condensates 
[ill l20ll2ll]. This effect is accounted for by the solution 
vs 7^ of Eq.© which becomes energetically favorable 
above Q CT . 

In the case of Bose-Einstein condensates the 
quadrupole instability occurs at f2 cr = u_i/\/2 when the 
lo = \/2uj± hydrodynamic quadrupole mode becomes en- 
ergetically unstable in the rotating frame [2(| • A different 
value is predicted to occur for the rotating Fermi gas, due 
to the new boundary conditions imposed by the presence 
of the normal component (22j. In fact, the current of the 
superfluid, in the rotating frame, should be tangential 
to the interface where the density, differently from the 
Bose-Einstein condensate (BEC) case, does not vanish. 

In the following we will determine the value of the crit- 
ical velocity in the simplest case of a rotating unpolarized 
gas and we will consider the onset of a quadrupole defor- 
mation. 

We consider an axially symmetric potential uj x = lo v = 
lu± =/= lo z and a solution where axi-symmetry is sponta- 
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neously broken. This is associated with the appearance 
of a non vanishing velocity field in the superfluid compo- 
nent whose velocity potential will be chosen of the form 

4> = ^/(r 2 ,z 2 ), (18) 

where / is an arbitrary function of r 2 = x 2 + y 2 and of 
z 2 . Note that in the case of the quadrupole instability of 
a Bose-Einstein condensate an exact stationary solution 
is found with / =const. 

The value of f2 cr corresponds to the onset of solutions 
with a deformed configuration. It is determined by solv- 
ing the continuity equation in the rotating frame 



V ■ [(V$ -fix r)n s (r)] = 0, 



(19) 



where ns(r) is now no longer axisymmctric and by im- 
posing that the superfluid current be tangential to the 
boundary B(r) of the superfluid 



(V$ - fi x r) ■ V.B(r) 



0. 



r=R s 



(20) 



The density, in the local density approximation, is 
given by 

» S M = ^(f) 3 VvW'*, (2D 



where V^ s (r) = V(r) 



mf!(r x vs)z is the ef- 



fective harmonic potential felt by the superfluid. By 
expanding the external potential to first order in / as 
V s (r) = V(r) + 5V s {r) with 5V s (r) = -mV$(n x r), 
one finds n s (r) = (fi- V{r) fl 2 - §(/x- V(r)) 1 / 2 -<5V s (r). 
At the same time the border can be written as B(r) = 
B°(r)+5B(r), where B°(r) is the radius of the superfluid 
given by Q 



R 2 (8) = ^ ( 



uj 2 , sin 2 i 
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1-7 



-fl 2 sin 2 



(22) 

and 8B(y) = — toV^ • (fi x r) is the linear perturbation 
due to the quadrupole symmetry breaking. 
Then, Eq 5.1(15]) and $kj§ reduce to 



(n-V(r))A<S>+^(V8V s (r))-(n x r)-^W(r)- V$ = 0, 

(23) 



and 



V$ • VB (r) - (fi x r) ■ VSB(r) 



--Rs(e) 



0. (24) 



By inserting the explicit expressions of the respective 
functions in Eqs. (f2"3]) and |24]) . after some straightfor- 
ward algebra [231 ] we obtain 



and 



(1 - 7 (1 - n 2 ) - 2f> 2 )/ + (1 - 7 (1 - fl 2 ))r 2 f r 



+ (l- 1 )z 2 f z (1-r 



2)3/2 



r=Rs(e) 



0, (26) 



where the latter equation is evaluated at the interface 
Eq.©. In Eqs.lgni H| / r(z) is the first derivative of / 
with respect to r (z), f IT r zz ) the second, and the expres- 
sions have been renormalized with respect to the radius 
of a superfluid at rest (i?g) 2 = 2/j,/muj 2 _ and we have 
made the substitution £l/u;± — > O and the assumption 
uj z = 

In the case of a two-dimensional system (Rs(0) = i?s) 
the previous equations can be easily solved. In this case 
indeed Eqs. (f25]) and ((26]) reduce to 



(2n 2 - l)/(r 2 ) + (2 - 3r 2 )/ r (r 2 ) + |(1 - r 2 )r 2 / rr (r 2 ) = 0, 

(27) 



and 



[l- 7 (l-^ 2 )-2fi 2 ]/(r 2 ) 
[l- 7 (l-^ 2 )]r 2 / r (r 2 ) 



r=R s 



the latter being evaluated at 



(1-7) 



[1-7(1-^ 2 )]' 



(28) 



(29) 



The solutions of the continuity equation Eq. (|2"T|) are 
hypergeometric functions 2F\{a, b, c; r 2 ) [24]], with the co- 
efficients a, b and c given by 

7 1 

a = - - -V25 + 480 2 , 
4 4 



y 2 

6 = - + -V25 + 480 2 , 
4 4 



(30) 



Inserting these solutions in the boundary condition 
((28]) . we find the value f2 cr = QA5uj± for the emergence 
of a spontaneous quadrupole deformation in two dimen- 
sions. 

For the three-dimensional case the calculation is more 
cumbersome and can be solved only numerically The 
results of the numerical calculation yield the estimate 
il cr ~ 0.5cjj^ for the critical angular velocity in three di- 
mensions. Notice that both in two and three dimensions 
the critical velocity is predicted to be smaller than the 
value Q cr = lj±/v / 2 holding in the BEC case. 



V. CONCLUSIONS 



(2fi 2 -l)/+|(l-r 2 -z 2 )(r 2 / rr + z 2 / zz ) 

+ (2 - 3r 2 - 2z 2 )/ r + i(l - r 2 - 4z 2 )f z = 0, (25) 



In conclusion we have analyzed the effect of adiabatic 
rotation on a polarized Fermi gas at unitarity assum- 
ing phase separation between a superfluid and a normal 
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phase. We find that the normal phase is energetically 
favoured by the rotation and thus the superfluid is fur- 
ther depleted with respect to the nonrotating configura- 
tion. The normal region exhibits the typical bulge ef- 
fect due to the centrifugal force while the superfluid is 
squeezed. This has clear observable effects on the den- 
sity profiles which can be addressed experimentally. A 
striking feature is that although the global polarization 
is not affected by the rotation, the concentration n.|/n-f 
at the border increases from the non-rotating value on 
the z axis to a maximum value in the xy plane. 

We have also addressed the question of quadrupole in- 
stability of the superfluid core, which produces a sponta- 
neous breaking of axial symmetry of the cloud. The crit- 
ical frequency for the onset of the instability turns out 
to be smaller than in the BEC case. Its measurement 
would provide a further crucial test of the mechanism 
of phase separation and of the equation of state of the 
normal phase [y, [j], l25[ . 

In our work we assume that the polarized system phase 
separates in only two phases. This assumption works ex- 
tremely well for the experiment carried out so far Q. 
In the rotating case other phases could show up and it 
would be very interesting to see how they affect our re- 
sults. For example, within BCS theory a third superfluid 
phase is found to occupy a small region at the interface 
[lOj]. Unfortunately we know that BCS theory is quan- 
titatively (and sometimes also qualitatively) not correct. 
A more microscopic investigation of the phase separation 
at unitarity would be necccssary to settle the problem. 

Finally, let us briefly argue on the possibility of hav- 
ing a rotating system without vortices in the superfluid, 



which is the main assumption of the present work. If the 
vortices could enter the superfluid, the lowest energy con- 
figuration would not be the one discussed in the paper, 
but rather a superfluid core with vortices sorrounded by 
a rotating normal shell, as also experimentally seen in |l[ . 
We expect that there exists a barrier for a vortex to en- 
ter the superfluid, as happens for a BE condensate. The 
interface caused by the polarized fcrmions could however 
change the scenario. On the one hand, the phase sepa- 
ration in the trap should disfavor the vortex formation 
as the superfluid density is finite at the interface. On 
the other hand, since there is a relative velocity between 
the superfluid and the normal shell, vortex nucleation 
could be favoured by a Kclvin-Hclmoltz-like mechanism. 
Moreover, we have shown that the interface reduces the 
critical angular velocity for a quadrupole instability. In 
the BEC case the latter is considered as a route toward 
vortex formation [ill ]. In the end, more theoretical and 
experimental work is needed to enlighten the issue of vor- 
tex nucleation in these polarized systems. 
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